Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales).

PREMISE: Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome. METHODS: We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω). RESULTS: Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK. CONCLUSIONS: We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.

Phylogenomics, plastome degradation and mycoheterotrophy evolution of Neottieae (Orchidaceae), with emphasis on the systematic position and Loess Plateau-Changbai Mountains disjunction of Diplandrorchis.

BACKGROUND: Mycoheterotrophy is a unique survival strategy adapted to dense forests and has attracted biologists' attention for centuries. However, its evolutionary origin and related plastome degradation are poorly understood. The tribe Neottieae contains various nutrition types, i.e., autotrophy, mixotrophy, and mycoheterotrophy. Here, we present a comprehensive phylogenetic analysis of the tribe based on plastome and nuclear ITS data. We inferred the evolutionary shift of nutrition types, constructed the patterns of plastome degradation, and estimated divergence times and ancestral ranges. We also used an integration of molecular dating and ecological niche modeling methods to investigate the disjunction between the Loess Plateau and Changbai Mountains in Diplandrorchis, a mycoheterotrophic genus endemic to China that was included in a molecular phylogenetic study for the first time. RESULTS: Diplandrorchis was imbedded within Neottia and formed a clade with four mycoheterotrophic species. Autotrophy is the ancestral state in Neottieae, mixotrophy independently originated at least five times, and three shifts from mixotrophy to mycoheterotrophy independently occurred. The five mixotrophic lineages possess all plastid genes or lost partial/all ndh genes, whereas each of the three mycoheterotroph lineages has a highly reduced plastome: one lost part of its ndh genes and a few photosynthesis-related genes, and the other two lost almost all ndh, photosynthesis-related, rpo, and atp genes. These three mycoheterotrophic lineages originated at about 26.40 Ma, 25.84 Ma, and 9.22 Ma, respectively. Diplandrorchis had presumably a wide range in the Pliocene and migrated southward in the Pleistocene. CONCLUSIONS: The Pleistocene climatic fluctuations and the resultant migration resulted in the Loess Plateau-Changbai Mountains disjunction of Diplandrorchis. In the evolution of mycoheterotrophic lineages, the loss of plastid-encoded genes and plastome degradation are staged and irreversible, constraining mycoheterotrophs to inhabit understories with low light levels. Accordingly, the rise of local forests might have promoted the origin of conditions in which mycoheterotrophy is advantageous.

The extremely reduced, diverged and reconfigured plastomes of the largest mycoheterotrophic orchid lineage.

BACKGROUND: Plastomes of heterotrophic plants have been greatly altered in structure and gene content, owing to the relaxation of selection on photosynthesis-related genes. The orchid tribe Gastrodieae is the largest and probably the oldest mycoheterotrophic clade of the extant family Orchidaceae. To characterize plastome evolution across members of this key important mycoheterotrophic lineage, we sequenced and analyzed the plastomes of eleven Gastrodieae members, including representative species of two genera, as well as members of the sister group Nervilieae. RESULTS: The plastomes of Gastrodieae members contain 20 protein-coding, four rRNA and five tRNA genes. Evolutionary analysis indicated that all rrn genes were transferred laterally and together, forming an rrn block in the plastomes of Gastrodieae. The plastome GC content of Gastrodia species ranged from 23.10% (G. flexistyla) to 25.79% (G. javanica). The plastome of Didymoplexis pallens contains two copies each of ycf1 and ycf2. The synonymous and nonsynonymous substitution rates were very high in the plastomes of Gastrodieae among mycoheterotrophic species in Orchidaceae and varied between genes. CONCLUSIONS: The plastomes of Gastrodieae are greatly reduced and characterized by low GC content, rrn block formation, lineage-specific reconfiguration and gene content, which might be positively selected. Overall, the plastomes of Gastrodieae not only serve as an excellent model for illustrating the evolution of plastomes but also provide new insights into plastome evolution in parasitic plants.

A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants.

Mycoheterotrophic and parasitic plants are heterotrophic and parasitize on fungi and plants, respectively, to obtain nutrients. Large-scale comparative genomics analysis has not been conducted in mycoheterotrophic or parasitic plants or between these two groups of parasites. We assembled a chromosome-level genome of the fully mycoheterotrophic plant Gastrodia elata (Orchidaceae) and performed comparative genomic analyses on the genomes of G. elata and four orchids (initial mycoheterotrophs), three parasitic plants (Cuscuta australis, Striga asiatica, and Sapria himalayana), and 36 autotrophs from various angiosperm lineages. It was found that while in the hemiparasite S. asiatica and initial mycoheterotrophic orchids, approximately 4-5% of the conserved orthogroups were lost, the fully heterotrophic G. elata and C. australis both lost approximately 10% of the conserved orthogroups, indicating that increased heterotrophy is positively associated with gene loss. Importantly, many genes that are essential for autotrophs, including those involved in photosynthesis, the circadian clock, flowering time regulation, immunity, nutrient uptake, and root and leaf development, were convergently lost in both G. elata and C. australis. The high-quality genome of G. elata will facilitate future studies on the physiology, ecology, and evolution of mycoheterotrophic plants, and our findings highlight the critical role of gene loss in the evolution of plants with heterotrophic lifestyles.

Comparative genomic analysis reveals metabolic flexibility of Woesearchaeota.

The archaeal phylum Woesearchaeota, within the DPANN superphylum, includes phylogenetically diverse microorganisms that inhabit various environments. Their biology is poorly understood due to the lack of cultured isolates. Here, we analyze datasets of Woesearchaeota 16S rRNA gene sequences and metagenome-assembled genomes to infer global distribution patterns, ecological preferences and metabolic capabilities. Phylogenomic analyses indicate that the phylum can be classified into ten subgroups, termed A-J. While a symbiotic lifestyle is predicted for most, some members of subgroup J might be host-independent. The genomes of several Woesearchaeota, including subgroup J, encode putative [FeFe] hydrogenases (known to be important for fermentation in other organisms), suggesting that these archaea might be anaerobic fermentative heterotrophs.

Nutrient-driven genome evolution revealed by comparative genomics of chrysomonad flagellates.

Phototrophic eukaryotes have evolved mainly by the primary or secondary uptake of photosynthetic organisms. A return to heterotrophy occurred multiple times in various protistan groups such as Chrysophyceae, despite the expected advantage of autotrophy. It is assumed that the evolutionary shift to mixotrophy and further to heterotrophy is triggered by a differential importance of nutrient and carbon limitation. We sequenced the genomes of 16 chrysophyte strains and compared them in terms of size, function, and sequence characteristics in relation to photo-, mixo- and heterotrophic nutrition. All strains were sequenced with Illumina and partly with PacBio. Heterotrophic taxa have reduced genomes and a higher GC content of up to 59% as compared to phototrophic taxa. Heterotrophs have a large pan genome, but a small core genome, indicating a differential specialization of the distinct lineages. The pan genome of mixotrophs and heterotrophs taken together but not the pan genome of the mixotrophs alone covers the complete functionality of the phototrophic strains indicating a random reduction of genes. The observed ploidy ranges from di- to tetraploidy and was found to be independent of taxonomy or trophic mode. Our results substantiate an evolution driven by nutrient and carbon limitation.

Integrative analysis of gene expression and alternative splicing in microalgae grown under heterotrophic condition.

Heterotrophic cultures are the most effective approach to overcome low growth rate challenge in the most commercial microalgae. However, the mechanism through which heterotrophic condition regulates algae metabolism are not completely clear. Alternative Splicing (AS) is a common posttranscriptional process by which transcriptome and proteome plasticity increases at different environmental conditions. To identify and characterize of AS events in Auxenochlorella protothecoides microalga grown in autotrophic and heterotrophic, RNA-Seq data were analysed. We found that AS increased with the transition from autotrophic to heterotrophic condition. 705 and 660 differentially expressed (DEG) and spliced (DAS) genes were identified for A.protothecoides was transferred from autotrophic to heterotrophic condition, respectively. Moreover, there was slight coverage between DEG and DAS genes. Furthermore, functional analysis showed that the DAS genes are most frequently related to ion binding and stimulus response. The results also indicated that prevalence of Intron retention is associated with down-regulation of the genes involved in carotenoid biosynthesis. This study provides valuable insights into transcriptional and posttranscriptional plasticity of microalgae during growth mode change.

Systems biology approach identifies functional modules and regulatory hubs related to secondary metabolites accumulation after transition from autotrophic to heterotrophic growth condition in microalgae.

Heterotrophic growth mode is among the most promising strategies put forth to overcome the low biomass and secondary metabolites productivity challenge. To shedding light on the underlying molecular mechanisms, transcriptome meta-analysis was integrated with weighted gene co-expression network analysis (WGCNA), connectivity analysis, functional enrichment, and hubs identification. Meta-analysis and Functional enrichment analysis demonstrated that most of the biological processes are up-regulated at heterotrophic growth condition, which leads to change of genetic architectures and phenotypic outcomes. WGNCA analysis of meta-genes also resulted four significant functional modules across logarithmic (LG), transition (TR), and production peak (PR) phases. The expression pattern and connectivity characteristics of the brown module as a non-preserved module vary across LG, TR, and PR phases. Functional analysis identified Carotenoid biosynthesis, Fatty acid metabolism and Methane metabolism as enriched pathways in the non-preserved module. Our integrated approach was applied here, identified some hubs, such as a serine hydroxymethyltransferase (SHMT1), which is the best candidate for development of metabolites accumulating strains in microalgae. Current study provided a new insight into underlying metabolite accumulation mechanisms and opens new avenue for the future applied studies in the microalgae field.

Unprecedented Parallel Photosynthetic Losses in a Heterotrophic Orchid Genus.

Heterotrophic plants are evolutionary experiments in genomic, morphological, and physiological change. Yet, genomic sampling gaps exist among independently derived heterotrophic lineages, leaving unanswered questions about the process of genome modification. Here, we have sequenced complete plastid genomes for all species of the leafless orchid genus Hexalectris, including multiple individuals for most, and leafy relatives Basiphyllaea and Bletia. Our objectives are to determine the number of independent losses of photosynthesis and to test hypotheses on the process of genome degradation as a result of relaxed selection. We demonstrate four to five independent losses of photosynthesis in Hexalectris based on degradation of the photosynthetic apparatus, with all but two species displaying evidence of losses, and variation in gene loss extending below the species level. Degradation in the atp complex is advanced in Hexalectris warnockii, whereas only minimal degradation (i.e., physical loss) has occurred among some "housekeeping" genes. We find genomic rearrangements, shifts in Inverted Repeat boundaries including complete loss in one accession of H. arizonica, and correlations among substitutional and genomic attributes. Our unprecedented finding of multiple, independent transitions to a fully mycoheterotrophic lifestyle in a single genus reveals that the number of such transitions among land plants is likely underestimated. This study underscores the importance of dense taxon sampling, which is highly informative for advancing models of genome evolution in heterotrophs. Mycoheterotrophs such as Hexalectris provide forward-genetic opportunities to study the consequences of radical genome evolution beyond what is possible with mutational studies in model organisms alone.

Genomic reconstructions and potential metabolic strategies of generalist and specialist heterotrophic bacteria associated with an estuary Synechococcus culture.

Interactions between photoautotrophs and heterotrophs are central to marine microbial ecosystems. Synechococcus are dominant marine phototrophs, and they are frequently associated with heterotrophic bacteria. These co-cultures provide a useful research system to investigate photoautotroph-heterotroph interactions in marine systems. Bacteria within the Roseobacter clade and Flavobacteria are two of the main bacterial lineages that exhibit intimate associations with Synechococcus populations. We conducted metagenomic analyses of a Synechococcus culture, followed by genomic binning of metagenomic contigs, and recovered five nearly complete genomes, including members of the Roseobacter clade (i.e. Marivita sp. XM-24) and Flavobacteria (i.e. Fluviicola sp. XM-24). Marivita sp. XM-24 is an ecological generalist of the Roseobacter clade and displays diverse metabolic capacities for the acquisition of nutrients and energy sources. Specifically, the genome contained numerous gene complements involved in the uptake and metabolism of nitrogen- and phosphorus-containing inorganic and organic compounds, in addition to the potential for aerobic anoxygenic photosynthesis, oxidation of carbon monoxide, inorganic sulfur oxidation, DMSP demethylation and PHA metabolism. The genome of the Flavobacteria representative, Fluviicola sp. XM-24, contained numerous peptidases, glycoside hydrolases, adhesion-related proteins and genes involved in gliding motility. Fluviicola sp. XM-24 likely specialize in the degradation of high molecular weight compound exudates from Synechococcus cells, including polysaccharides and polypeptides via attachment to particles, surfaces or cells. The distinct metabolic strategies identified within several heterotrophic bacteria that are associated with Syneochococcus cells provide insights into their lifestyles and nutrient utilization patterns, in addition to their interactions with photoautotrophs. Biological interactions, including mutualism, competition and antagonism, shape the microbial community structure of marine environments and are critical for understanding biogeochemical cycling in the ocean. These results provide valuable insights into the nature of interactions between dominant marine photoautotrophs and associated bacterial heterotrophs.

Effects of carbon source, C/N ratio, nitrate, temperature, and pH on N2O emission and functional denitrifying genes during heterotrophic denitrification.

The effects of operational parameters such as carbon source, C/N ratio, initial nitrate concentration, temperature, and pH value on heterotrophic denitrification and functional denitrifying genes were evaluated. When methanol was used as the sole carbon source, complete denitrification was performed in a short time without nitrous oxide (N2O) emission. Complete denitrification was performed at high C/N ratios (5.14 and 12.85) and low initial nitrate concentrations (75.9 and 151.6 mg N L-1). The denitrification rate was not temperature-sensitive in the range of 25-35 °C, but tended to decrease at a low pH of 5-6. The relationships between N2O emission and functional genes under various operational conditions were investigated by Pearson correlation and association network analyses. The C/N ratio was a key factor for N2O emission during the heterotrophic denitrification process. This information on the denitrification performance and its association with functional gene dynamics under various operational conditions is useful for N2O mitigation strategies for wastewater treatment processes.

Optimization of carbon and energy utilization through differential translational efficiency.

Control of translation is vital to all species. Here we employ a multi-omics approach to decipher condition-dependent translational regulation in the model acetogen Clostridium ljungdahlii. Integration of data from cells grown autotrophically or heterotrophically revealed that pathways critical to carbon and energy metabolism are under strong translational regulation. Major pathways involved in carbon and energy metabolism are not only differentially transcribed and translated, but their translational efficiencies are differentially elevated in response to resource availability under different growth conditions. We show that translational efficiency is not static and that it changes dynamically in response to mRNA expression levels. mRNAs harboring optimized 5'-untranslated region and coding region features, have higher translational efficiencies and are significantly enriched in genes encoding carbon and energy metabolism. In contrast, mRNAs enriched in housekeeping functions harbor sub-optimal features and have lower translational efficiencies. We propose that regulation of translational efficiency is crucial for effectively controlling resource allocation in energy-deprived microorganisms.

Trait-based patterns of microbial dynamics in dormancy potential and heterotrophic strategy: case studies of resource-based and post-press succession.

Understanding the relationship between microbial community structure and function is a major challenge in microbial ecology. Recent work has shown that community weighted mean 16S rRNA gene copies, as a proxy for heterotrophic growth strategy, is a microbial community trait that decreases predictably over successional trajectories that are underpinned by changes in resource availability. However, it has been challenging to identify other microbial traits that are predictive of community functions and have consistent patterns with succession. Trait-based patterns of secondary succession (e.g., after a disturbance) are less often considered, and these responses may be underpinned by abiotic drivers other than changes in resources. In this perspectives piece, we present hypotheses about microbial traits important for microbial succession in resource-based and post-press disturbance scenarios, as synthesized from previous works and extended within this work. Using four case studies, we compare two traits, heterotrophic strategy and dormancy potential, and two different types of succession, resource-based (endogenous heterotrophic) and post-press. There were decreases in weighted ribosomal operon counts and in dormancy genes over resource-based succession. Both traits also were lower in post-press succession as compared to reference conditions, but increased with time from disturbance. Thus, dormancy potential may be an additional trait that changes predictably with succession. Finally, considering changes in microbial community traits over post-press succession is as important as over resource-based succession. These patterns need to be interpreted carefully and reference and recovering samples can be collected to improve interpretation of changes in community traits over post-press succession.

Phylogenomic inference in extremis: A case study with mycoheterotroph plastomes.

PREMISE OF THE STUDY: Phylogenomic studies employing large numbers of genes, including those based on plastid genomes (plastomes), are becoming common. Nonphotosynthetic plants such as mycoheterotrophs (which rely on root-associated fungi for essential nutrients, including carbon) tend to have highly elevated rates of plastome evolution, substantial genome reduction, or both. Mycoheterotroph plastomes therefore provide excellent test cases for investigating how extreme conditions impact phylogenomic inference. METHODS: We used parsimony and likelihood analysis of protein-coding gene sets from published and newly completed plastomes to infer the phylogenetic placement of taxa from the 10 angiosperm families in which mycoheterotrophy evolved. KEY RESULTS: Despite multiple very long branches that reflect elevated substitution rates, and frequently patchy gene recovery due to genome reduction, inferred phylogenetic placements of most mycoheterotrophic lineages in DNA-based likelihood analyses are both well supported and congruent with other studies. Amino-acid-based likelihood placements are broadly consistent with DNA-based inferences, but extremely rate-elevated taxa can have unexpected placements-albeit with weak support. In contrast, parsimony analysis is strongly misled by long-branch attraction among many distantly related mycoheterotrophic monocots. CONCLUSIONS: Mycoheterotrophic plastomes provide challenging cases for phylogenomic inference, as substitutional rates can be elevated and genome reduction can lead to sparse gene recovery. Nonetheless, diverse likelihood frameworks provide generally well-supported and mutually concordant phylogenetic placements of mycoheterotrophs, consistent with recent phylogenetic studies and angiosperm-wide classifications. Previous predictions of parallel photosynthesis loss within families are supported for Burmanniaceae, Ericaceae, Gentianaceae, and Orchidaceae. Burmanniaceae and Thismiaceae should not be combined as a single family in Dioscoreales.

Genome Reports: Contracted Genes and Dwarfed Plastome in Mycoheterotrophic Sciaphila thaidanica (Triuridaceae, Pandanales).

With a reduced need for photosynthesis, the plastome of parasitic and mycoheterotrophic plants degrades. In the tiny, fully mycoheterotrophic plant Sciaphila thaidanica, we find one of the smallest plastomes yet encountered. Its size is just 12,780 bp and it contains only 20 potentially functional housekeeping genes. Thus S. thaidanica fits the proposed model of gene loss in achlorophyllous plants. The most astonishing feature of the plastome is its extremely compact nature, with more than half of the genes having overlapping reading frames. Additionally, intergenic sequences have been reduced to a bare minimum, and the retained genes have been reduced in length both compared with the orthologous genes in another mycoheterotrophic species of Sciaphila and in the autotrophic relative Carludovica.

Evolution of heterotrophy in chrysophytes as reflected by comparative transcriptomics.

Shifts in the nutritional mode between phototrophy, mixotrophy and heterotrophy are a widespread phenomenon in the evolution of eukaryotic diversity. The transition between nutritional modes is particularly pronounced in chrysophytes and occurred independently several times through parallel evolution. Thus, chrysophytes provide a unique opportunity for studying the molecular basis of nutritional diversification and of the accompanying pathway reduction and degradation of plastid structures. In order to analyze the succession in switching the nutritional mode from mixotrophy to heterotrophy, we compared the transcriptome of the mixotrophic Poterioochromonas malhamensis with the transcriptomes of three obligate heterotrophic species of Ochromonadales. We used the transcriptome of P. malhamensis as a reference for plastid reduction in the heterotrophic taxa. The analyzed heterotrophic taxa were in different stages of plastid reduction. We investigated the reduction of several photosynthesis related pathways e.g. the xanthophyll cycle, the mevalonate pathway, the shikimate pathway and the tryptophan biosynthesis as well as the reduction of plastid structures and postulate a presumable succession of pathway reduction and degradation of accompanying structures.

Abundance and diversity of aerobic heterotrophic microorganisms and their interaction with cyanobacteria in the oxic layer of an intertidal hypersaline cyanobacterial mat.

Aerobic heterotrophic microorganisms (AH) play a significant role in carbon cycling in cyanobacterial mats; however, little is known about their abundance, diversity and interaction with cyanobacteria. Using catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH), bacterial counts in the mat's oxic layer reached a mean of 2.23 ± 0.4 × 1010 cells g-1. Cultivation of AH yielded strains belonging to Actinobacteria, Bacteroidetes, Firmicutes, Gammaproteobacteria and Haloarchaea. 16S rRNA bacterial sequences retrieved from the mat's oxic layer were related to Bacteroidetes, Chloroflexi and Proteobacteria, whereas archaeal sequences belonged to Crenarchaeota and Haloarchaea. Monocultures of cyanobacteria from the same mat were associated with different AH, although Bacteroidetes were found in most cultures. CARD-FISH showed that Bacteroidetes- and Chloroflexi-related bacteria were closely associated with filaments of Microcoleus chthonoplastes. The growth of an axenic culture of M. chthonoplastes PCC7420 was stimulated on the addition of a filtrate obtained from a non-axenic Microcoleus culture and containing only AH and released substances. In contrast, a similar filtrate from a non-axenic Cyanothece-related culture killed Cyanothece PCC 7418. We conclude that a diverse community of AH exist in close association with cyanobacteria in microbial mats and the interactions between AH and cyanobacteria are species-specific and involve the release of substances.

Activity and functional properties of the isocitrate lyase in the cyanobacterium Cyanothece sp. PCC 7424.

Cyanobacteria are ubiquitous photoautotrophs that assimilate atmospheric CO2 as their main source of carbon. Several cyanobacteria are known to be facultative heterotrophs that are able to grow on diverse carbon sources. For selected strains, assimilation of organic acids and mixotrophic growth on acetate has been reported for decades. However, evidence for the existence of a functional glyoxylate shunt in cyanobacteria has long been contradictory and unclear. Genes coding for isocitrate lyase (ICL) and malate synthase were recently identified in two strains of the genus Cyanothece, and the existence of the complete glyoxylate shunt was verified in a strain of Chlorogloeopsis fritschii. Here, we report that the gene PCC7424_4054 of the strain Cyanothece sp. PCC 7424 encodes an enzymatically active protein that catalyses the reaction of ICL, an enzyme that is specific for the glyoxylate shunt. We demonstrate that ICL activity is induced under alternating day/night cycles and acetate-supplemented cultures exhibit enhanced growth. In contrast, growth under constant light did not result in any detectable ICL activity or enhanced growth of acetate-supplemented cultures. Furthermore, our results indicate that, despite the presence of a glyoxylate shunt, acetate does not support continued heterotrophic growth and cell proliferation. The functional validation of the ICL is supplemented with a bioinformatics analysis of enzymes that co-occur with the glyoxylate shunt. We hypothesize that the glyoxylate shunt in Cyanothece sp. PCC 7424, and possibly other nitrogen-fixing cyanobacteria, is an adaptation to a specific ecological niche and supports assimilation of nitrogen or organic compounds during the night phase.

Plastomes on the edge: the evolutionary breakdown of mycoheterotroph plastid genomes.

Contents 48 I. 48 II. 50 III. 53 54 References 54 SUMMARY: We examine recent evidence for ratchet-like genome degradation in mycoheterotrophs, plants that obtain nutrition from fungi. Initial loss of the NADH dehydrogenase-like (NDH) complex may often set off an irreversible evolutionary cascade of photosynthetic gene losses. Genes for plastid-encoded subunits of RNA polymerase and photosynthetic enzymes with secondary functions (Rubisco and ATP synthase) can persist initially, with nonsynchronous and quite broad windows in the relative timing of their loss. Delayed losses of five core nonbioenergetic genes (especially trnE and accD, which respectively code for glutamyl tRNA and a subunit of acetyl-CoA carboxylase) probably explain long-term persistence of heterotrophic plastomes. The observed range of changes of mycoheterotroph plastomes is similar to that of holoparasites, although greater diversity of both probably remains to be discovered. These patterns of gene loss/retention can inform research programs on plastome function.

Lineage-Specific Reductions of Plastid Genomes in an Orchid Tribe with Partially and Fully Mycoheterotrophic Species.

The plastid genome (plastome) of heterotrophic plants like mycoheterotrophs and parasites shows massive gene losses in consequence to the relaxation of functional constraints on photosynthesis. To understand the patterns of this convergent plastome reduction syndrome in heterotrophic plants, we studied 12 closely related orchids of three different lifeforms from the tribe Neottieae (Orchidaceae). We employ a comparative genomics approach to examine structural and selectional changes in plastomes within Neottieae. Both leafy and leafless heterotrophic species have functionally reduced plastid genome. Our analyses show that genes for the NAD(P)H dehydrogenase complex, the photosystems, and the RNA polymerase have been lost functionally multiple times independently. The physical reduction proceeds in a highly lineage-specific manner, accompanied by structural reconfigurations such as inversions or modifications of the large inverted repeats. Despite significant but minor selectional changes, all retained genes continue to evolve under purifying selection. All leafless Neottia species, including both visibly green and nongreen members, are fully mycoheterotrophic, likely evolved from leafy and partially mycoheterotrophic species. The plastomes of Neottieae span many stages of plastome degradation, including the longest plastome of a mycoheterotroph, providing invaluable insights into the mechanisms of plastome evolution along the transition from autotrophy to full mycoheterotrophy.